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ABSTRACT 

We report the discovery of a transiting planet orbiting the star TYC 6446-326-1. The 
star, WASP-22, is a moderately bright (V=12.0) solar-type star (Teff 6000 ± lOOK, [Fe/H] = 
—0.05 ± 0.08). The lightcurve of the star obtained with the WASP-South instrument shows pe- 
riodic transit-like features with a depth of about 1% and a duration of 0.14 d. The presence of a 
transit-like feature in the lightcurve is confirmed using z-band photometry obtained with Faulkes 
Telescope South. High resolution spectroscopy obtained with the CORALIE and HARPS spec- 
trographs confirm the presence of a planetary mass companion with an orbital period of 3.533 days 
in a near-circular orbit. From a combined analysis of the spectroscopic and photometric data 
assuming that the star is a typical main-sequence star we estimate that the planet has a mass 
Mp — 0.56 ± 0.02Afjup and a radius Rp = 1.12 ± 0.04i?jup. In addition, there is a linear trend of 
40ms~^ in the radial velocities measured over 16 months, from which we infer the presence 
of a third body with a long period orbit in this system. The companion may be a low mass 
M-dwarf, a white dwarf or a second planet. 
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Introduction 



The WASP project (jPoUacco et al.ll2006[ ) is cur- 
rently one of the most successful wide-area surveys 
designed to find exoplanets transiting bright stars 
(V < 12.5). Other succe sful sur veys include HAT- 
net (iBakos et alj |2004|). XO dMcCuUough et al 



I2OO5I ) and TrES (|0'Donovan et al.ll2006l ). There is 
continued interest in finding transiting exoplanets 
because they can be accurately characterized and 
studied in some detail, e.g., the mass and radius of 
the planet can be accurately measured. This gives 
us the opportunity to explore the relationships be- 
tween the density of the planet and other proper- 
ties of the planetary system, e.g., the semi- major 
axis, the spectral type of the star, the eccentric- 
ity of the orbit, etc. Given the wide variety of 
transiting planets being discovered and the large 
number of parameters that characterize them, sta- 
tistical studies will require a large sample of sys- 
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terns to identify and quantify the relationships be- 
tween these parameters. These relationships can 
be used to test models of the formation, structure 
and evolution of short period exoplanets. 

A particular puzzle related to the properties 
of hot Jupiters is the wide range in their densi- 
ties. Very dense hot Jupiters such as HD 149026 
are though t to contain a dense, metallic core 



Table 1: Radial velocity measurements. 



(jSato et al.l 120051 ). There is currently no gen- 
erally agreed explanation for the existence of 
hot Jupiters with densities 5-10 times lower 
than the density of Jupiter, e .g. WASP-17b 
( Anderson et al.l2010[ ). TrES-4b (iMandushev et"al 
2007|) and WASP-12b (|Hebb et al.l l2009l ). One" 
possibility is that the planets are heated by tidal 
forces, and that these are driven by the presenc e 
of a third body in the system ( Mardlina 12007 ). 
Other possibiliti es include enhanced opacity in 
the atmosphere (jBurrows et al.l 120071) . the distri- 
butio n of heavy elements in the core ( jBaraffe et al 



20081 ) and kinetic heating f rom the irradiated at- 
mosp here into the interior ([Showman and Guillot 
20021 ). 

Here we report the discovery of a hot Jupiter 
system, WASP-22, identified using the WASP- 
South instrument and present evidence that it is 
a member of a hierarchical triple system. 

2. Observations 



Th e WA SP survey is describe d in lPoUacco et al 
(|2006l ) and lWilson et al.l (|2008l) while a discussion 
of our candidate selection methods can be found 



- . Collier Cameron et al.l (|2007l ). iPoUacco et al 
( 20081 ). and references therein. 



The WASP-South instrument consists of 8 
cameras, each with a Canon 200-mm f/1.8 lens 
and a 2kx2k e2V CCD detector resulting in 
an image scale of approximately 14 arcsec per 
pixel. The star TYC 6446-326-1 ISWASP 
J033116.32-234911.0) was observed 3133 times 
by one camera on the WASP-South instrument 
from 2006 August to 2007 January. A further 
6282 obervations were obtained with the same 
camera from 2007 August to 2008 January. The 
star also appeared in the images obtained with a 
second camera during the second observing sea- 
son, so a further 5889 observations were obtained 
with this camera during that interval. 

The WASP-South lightcurves of WASP-22 
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Fig. 1.— WASP-South photometry of WASP-22 
folded on the orbital period P = 3.532759 d. Up- 
per panel: all data. Lower panel: data within 0.12 
phase units of mid-transit together with the model 
fit described in Section [XT] (solid line). 
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Fig. 2. — Faulkes Telescope South z-band pho- 
tometry of WASP-22 (points) with the model fit 
described in Section [01 (solid line). 



show transit-like features with a depth of approx- 
imately 0.012 magnitudes recurring with a 3.53-d 
period (Fig. [T]) . These were detected in the WASP 
photometry from the 2006 season using the de- 
tr ending and transit detection m ethods described 



Collier Cameron et al.l (|2007l ). but judged to be 



"unconvincing" at that stage. The same 3.53-d pe- 
riodicity was detected in the 2007 data from both 
cameras, which was taken as strong evidence that 
the periodic transit signal was real. The spectral 
type of the star was estimated to be approximately 
Gl based on the catalogue photometry available 
for this star at the time. The duration and depth 
of the transit is consistent with the hypothesis 
that it is due to the transit of a planet-like com- 
panion to a main-sequence Gl star and the WASP 
lightcurves show no indication of any ellipsoidal 
variation due to the distortion of the star by a 
massive companion. 

We also considered whether WASP-22 was 
likely to be a dwarf star given its J— H colour 
and J-band reduced proper motion, Hj. The 
position of WASP-22 relative to the polyno- 
mi al boundary between d w arfs a nd giants given 

(|2007l) 



Collier Cameron et al. 



suggests that 
WASP-22 is a dwarf. We have also used the 
results of the RAVE survey DR2 (jZwitter et al 



2008[ ) to calibrate the probability that a star is 
a dwarf given the observed values of J— H and 
Hj. The RAVE survey observed stars at a similar 
range of magnitudes and galactic positions to the 
WASP survey and so it is well suited to this pur- 
pose. Stars listed in the survey with \ogg > 3.5 
were identified as dwarf stars, those with lower 
\ogg estimates were identified as giants. We then 
created a look-up table of the relative numbers of 
dwarfs and giants as a function of position in the 
Hj v. J— H plane. This criterion also suggested 
that WASP-22 is likely to be a dwarf star because 
the ratio of dwarfs to giants in this region of the 
Hj V. J— H plane is dwarf:giant = 10:0. 

The star was included as a "grade-A" candi- 
date in our programme of follow-up spectroscopic 
observations using the CORALIE spectrograph on 
the Euler 1.2-ni telescope and the HARPS spectro- 
graph on the 3.6-m ESO telescope, both located 
at La Silla, Chile. We obtained 37 radial-velocity 
measurements during the interval 2008 August 27 
to 2009 December 25 with CORALIE and 6 mea- 
surements with HARPS during the interval 2008 
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Oct 4-16. The CORALIE spectra have a signal- 
to-noise ratio of 10-20, while the typical signal- 
to-noise ratio of the HARPS spectra is 50. The 
measurements are given in Table [TJ where we also 
provide the bisector span, BS, which measures the 
asymmetry of the cross-correlation function. The 
standard error of the the bisector span measure- 
ments is 2(Trv- 

We also obtained photometry of TYC 6446-326- 
1 and other nearby stars on 2009 August 8 using 
the LCOGT 2.0-m Faulkes Telescope South (FTS) 
at Siding Spring Observatory. The Merope cam- 
era we used has an image scale of 0.279 arcsec- 
onds/pixel when used in the 2x2 binning mode we 
employed. We used a Pan-STARR^I z-band filter 
to obtain 107 images covering one ingress of the 
transit. These images were processed in the stan- 
dard way with irafI using a stacked bias image, 
dark frame, and sky flat. Minimal fringing was 
present in the z-band images due to the deep de- 
pletion CCD in the camera, so no fringe correction 
was applied. The DAOPHOT photometry pack- 
age (|Stetsonlll987l ) was used to perform object de- 
tection and aperture photometry with an aperture 
size of 10 binned pixels in radius. The 5' x 5' field 
of view of the instrument contained 6 comparison 
stars that were used in deriving the differential 
magnitudes with a photometric precision of 3.1 
mmag. The coverage of the out-of-transit phases 
is quite limited, but the data are sufficient to con- 
firm that transit-like features seen in the WASP- 
South data are due to the star TYC 6446-326-1 
and to provide better measurements of the depth 
of the transit and the duration of ingress than is 
possible from the WASP-South data (Fig.^]). We 
note that the star ISWASP J033121. 40-234857.8 
77arcsec from WASP-22 is also a variable star. 
This star showed a dip in brightness of about 
0.01 magnitudes lasting about 100 minutes during 
our observation of WASP-22. 

We examined the power spectra of the WASP 
data from each camera and season separately to 
look for any intrinsic photometric variability due 
to magnetic activity in WASP-22. There is no pe- 
riodic signal consistently detected in all three data 



^http: //pan- Starrs . if a.hawaii . edu/public" 

■^IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 



sets SO any such variability must have an ampli- 
tude less than a few milli-magnitudes. 

All photometric presented in this paper are 
available from the NStED databaseH 

3. WASP-22 Stellar Parameters 

The individual HARPS spectra of WASP-22 
were co-added to produce a single spectrum with 
a typical S/N of around 100:1 that we have anal- 
ysed to determine the atmospheric parameters of 
the star. The standard pipeline reduction prod- 
ucts were used in the analysis. 

The analysis was performed using t he UCLSYN 
spect ral synthesis package (Smith 1992; jSmallev et al 
20011 ) and ATLASQ models with out convective over- 
shooting ( Castelli et al.|[l997h . The Ha line was 
used to determine the effective temperature (Toff), 
while the Na i D and Mg I b lines were used as sur- 
face gravity (logt;) diagnostics. The parameters 
obtained from the analysis are listed in Table [5J 

The equivalent widths of several clean and un- 
blended lines were meas ured. Atomic line data 
was mainly taken from the lKurucz and "Beiil (|l995h 
compilation, but with updated van der Waals 
broadening coefficients for lin es in ^Barkle m et al 



2000 ) a nd log qf values f romlGonzalez and Laws 



2000 ). iGonzalez et"al] (|200l[ ) or ISantos et al 



20041 ). A value for microturbulence (^t) was de- 



termin ed from Fe i using the method of iMagainl 
(|l984l ). The ionization balance between Fe I and 
Fe II and the null-dependence of abundance on ex- 
citation potential wer e used as an a dditional Toff 
and logg diagnostics ( Smallev 20051 ). 

We have determined the elemental abundances 
of several elements (also listed in Table ^ from 
their measured equivalent widths. The quoted er- 
ror estimates include that given by the uncertain- 
ties in Tcff, log 5 and ^t, as well as the scatter 
due to measurement and atomic data uncertain- 
ties. There is no evidence from these data for any 
abundance anomalies in WASP-22. 

The spectrum of WASP-22 shows a clear Li i 
6708A line with EW = 32 ± 1 mA, indicating 
an abundance of log n(Li/H) -f 12 = 2.23 ± 0.08 
dex. We have compared this lithium abundance 
to the relations between effective t emperature and 
age bv lSestito and Randich ( 2005 ). For stars with 



• ^http : //nsted ■ Ipac . caltech ■ edu | 
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Fig. 3. — Radial velocity and bisector span mea- 
surements for WASP-22. Upper panel: Radial ve- 
locity relative to the centre-of-mass velocity in- 
cluding the long-term drift ^ (points with error 
bars) compared to our model for the spectroscopic 
orbit (solid line). Lower panel: bisector span mea- 
surements. Data obtained with the HARPS spec- 
trograph are plotted with filled circles, CORALIE 
data with open circles. 
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Fig. 4. — Residuals from a Keplerian orbit fit to 
our radial velocity data without a linear trend in 
the model. Symbols are as for Fig. [3l 



Tefi « 6000 K we find that this lithium abundance 
implies a lower limit to the age of about 1 Gyr. 

The projected stellar rotation velocity (wsini) 
was determined by fitting the profiles of several 
unblended Fe i lines. A value for macroturbulence 
(wmac) of 4.5 km s'~^ was assumed, based on the 
tabulation by iGravl (|2008h . and an instrumental 
FWHM of 0.065A, determined from the telluric 
lines around 6300A. A best fitting value of wsini 
= 3.5 ± 0.6 kms~^ was obtained. Inspection of 
the Call H and K lines shows no indication of 
any emission due to chromospheric activity, as ex- 
pected for such a slowly rotating star. If we as- 
sume that the rotation axis of the star is approxi- 
mately aligned with the orbital axis of the planet's 
orbit we can estimate that the rotation period of 
the star is fin t = 16 ± 3 d. The rota tion-colour- 
age relation of ICollier Cameron et al.l ([2009) then 
yields an age estimate of approximately 3 ± 1 Gyr 
for this star. We have also compared the values 
of Teff and th e stella r density p^, to the models of 
Girardi et al. (l2000l) and find that the age of the 
star is likely to be less than 6 Gyr based on these 
models. 

We did not see any indication of additional 
spectral lines in the spectrum nor any trend in 
equivalent widths that might suggest contamina- 
tion of the spectrum by another star. We estimate 
that the third body in this system discussed below 
contributes less than 5% of the light in the optical 
spectrum. 

In summary, WASP-22 appears to be a slowly 
rotating, chromospherically inactive, main-sequence 
star that is slightly hotter than the Sun. 

3.1. Planetary parameters 

The amplitude of the radial velocity variation 
with the same period as the transit lightcurve 
(Fig. [3]) and the lack of any correlation between 
this variation and the bisector span establish the 
prese nce of a planetary mass companion to this 
star (jQueloz et al.ll200l[ ). 

The CORALIE and HARPS radial velocity 
measurements were combined with the WASP- 
South and FTS photometry in a simultaneous 
Markov-chain Monte-Carlo (MCMC) analysis to 
find the parameters of the WASP-22 system. The 
shape of the transit is not well defined in the 
WASP-South or FTS photometry, so we have im- 
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posed an assumed main-sequence mass-radius re- 
lation as an additional constraint in our anal- 
ysis of the data. The stellar mass is deter- 
mined from the parameters Tgg, log o and [Fe/H] 



using the procedure described by lEnoch et al 
()2010l) . based on the compi lation of eclipsing 



binary data by Irorres et all (|2010[ ). The code 
uses Tcff and [Fe/H] as MCMC jump variables, 
constrained by Bayesian priors based on the 
spectroscopically-determined values given in Ta- 
blej2l Limb-darkening coefficients are taken from 
Claretl (UoOO). 



A careful analysis of the radial velocity data 
clearly shows a trend in the residuals as a func- 
tion of time (Fig. S]). We therefore adapted our 
MCMC analysis to include an extra parameter, 
which describes a linear trend in the centre- 
of-mass velocity of the star. 

The parameters derived from our MCMC anal- 
ysis are listed in Table |31 We found that we did 
not need to account for any additional noise in the 
radial velocity data to due to stellar activity ("jit- 
ter"). The contribution to the total chi-squared 
from our 43 radial velocity measurements is 32.4. 
The surface gravity derived from our MCMC so- 
lution is consistent with the log 5 value from the 
analysis of the spectrum, but the large uncertainty 
on the latter value means that this is a rather weak 
constraint. 

We considered the contribution of red noise to 
the standard errors quoted in Table [3] While 
the individual transits in the WASP data are af- 
fected by red noise, the analysis of the combined 
lightcurve covering many individual transits will 
not be strongly affected by red noise because there 
will be no correlation between the systematic noise 
from different nights. The FTS lightcurve is af- 
fected by red noise so we have investigated the 
effect of this using the "prayer bead" method. A 
separate MCMC analysis was performed in which 
synthetic FTS lightcurves were created from the 
model fit to the lightcurve and the residuals from 
this model after cyclic permutation at each step 
in the MCMC chain. We find that this does not 
significantly increase the error estimates for any 
of the parameters with the exception of the tran- 
sit epoch and period. The standard errors for 
these two parameters in Table |3] are taken from 
the "prayer bead" MCMC analysis. 

We have included the parameters e sin(aj) and 



e cos(w) as free parameters in the solution because 
the orbit of the planet is not known to be circular, 
although the values of derived and their standard 
errors are quite consistent with this hypothesis. 
Imposing a circular orbit in the solution has a neg- 
ligible effect on the value of the other parameters 
derived, but would reduce the estimated errors on 
these parameters, perhaps unrealisitically so. 

4. Discussion 

Apart from the presence of a third body, 
WASP-22 is a typical hot Jupiter planetary 
system. WASP-22 b joins a growing number 
of planets discovered with masses ~ 0.5A/jup, 
radii ~ li?jup orbiting solar-like stars with pe- 
riods of about 3 days, e.g., WASP-13b, WASP- 
llb/HAT-P-lOb, CoRoT-5b, WASP-6b, HAT- 
P-lb, OGLE-TR-lllb, WASP-25b, HAT-P-3b, 
Kepler-8 b, ctcEl 

The third body we have detected from the lin- 
ear trend in the radial velocities contributes less 
than 5% of the fiux in the optical spectrum, so 
it is at least 3 magnitudes fainter than WASP- 
22 at V. This rules out the possibilty that the 
companion is K-type dwarf star, but leaves the 
possibility that the companion is an M-dwarf or 
a white dwarf. There is good agreement be- 
tween the effective temperature we derive from 
the analysis of the optical spectrum and the ef- 
fective t emperature derived us ing the infrared flux 



method ( Blackwell et al. I I1979I) . We estimate that 



we would have been able detect a cool compan- 
ion to WASP-22 if the 2MA SS K.,-band magni - 



tude (K, = 10.318±0.020. ISkrutskie et al.ll2006t ) 

were increased by more than 0.1 magnitudes. This 
also rules out a K-dwarf companion but not an 
M-dwarf companion or a white dwarf companion. 
For example, a K-dwarf companion to WASP-22 
would produce K-band infrared excess of 0.1 mag- 
nitudes or more, in which case the effective tem- 
perature derived using the infrared flux method 
would be inconsistent with the value derived from 
the analysis of the optical spectrum. 

There are few direct constraints on the prop- 
erties of the third body from the radial velocity 
data available to date. The observation that the 
trend is linear over 16 months suggests that the 



' ^http://exoplanet.eu\ 
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orbital period is at least a few years. If the or- 
bit of the third body is approximately circular 
and coplanar with the inner orbit with a period 
of several years and an amplitude comparable to 
the velocity range observed so far, then it is pos- 
sible that the third body is a second planet, i.e., 
a configuration similar to the do uble-planet sys- 
tem HAT-P-13 (|Bakos et al.ll2009l ). Further radial 
velocity monitoring will be required to determine 
whether any or all of these assumptions is reason- 
able. This is certainly worth doing because the de- 
tection of a second planet in the WASP-22 system 
may make it possible to infer the int erior structure 



of WA SP-22 b using the method of iBatvgin et al 
()2009l) . This method relies on an accurate mea- 



surement of the eccentricity for the orbit of the 
inner planet. The data available to date are only 
of sufficient quality to state that the eccentricity 
of the orbit is small (e < 0.06) so more data of the 
quality we obtain from HARPS will be required 
to measure a precise value for the eccentricity of 
WASP-22 b's orbit. 

The distance to WASP-22 is approximately 
300pc so it may be possible to detect the compan- 
ion directly if it is a M-dwarf using high-constrast, 
high-resolution imaging. 

5. Conclusion 

The star WASP-22 (TYC 6446-326-1) has a hot 
Jupiter companion. A long-term linear trend in 
the mean value of the radial velocity shows that 
WASP-22 has a distant companion, i.e., it is a hier- 
archical triple system. The properties of the third 
body are poorly constrained by the data available 
to-date, but it may be an M-dwarf, a white dwarf 
or a second planet. 

WASP-South is hosted by the South African 
Astronomical Observatory and we are grateful for 
their ongoing support and assistance. Funding 
for WASP comes from consortium universities and 
from the UKs Science and Technology Facilities 
Council. 
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Table 3: System parameters for WASP-22. 
value for the Bond albedo A= 0. N.B. an 
ad ditional const r aint i n this solution so the 
see lEnoch et all (|2010t)') for details. 



The planet equilibrium temperature is calculated assuming a 
assumed main-sequence mass-radius relation is imposed as an 
mass and radius of the star are not independent parameters - 



Parameter 



Symbol 



Value 



Unit 



Transit epoch (BJD) 

Orbital period 

Transit duration 

Planet/star area ratio 

Impact parameter 

Stellar reflex velocity 

Centre-of-mass velocity at time Tq 

Velocity offset, HARPS - CORALIE 

Drift in centre-of-mass velocity 

Orbital separation 

Orbital inclination 

Orbital eccentricity 

Arg. of periastron 



Stellar mass 
Stellar radius 
Stellar surface gravity 
Stellar density 
Planet mass 
Planet radius 
Planet surface gravity 
Planet density 
Planet equil. temp. 



P 

Tu 

Rp/Rl 
b 

Ki 

7coralie(Jc) 
7harps ~ 7coralie 

It 

a 
i 
e 

ecos(w) 
e sin(ti;) 

R* 

log g* 

P* 

Mp 

Rv 

loggp 

pp 
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70.0 ± 1.7 
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40 ±5 
0.0468 ± 0.0004 
89.2 ±0.5 
0.023 ±0.012 

97-1-51 
^'-78 

0.012 ±0.011 
0.006 ±0.021 

1.1 ±0.3 
1.13 ±0.03 
4.37 ±0.02 

0.76±0.06 
0.56 ±0.02 
1.12 ±0.04 
3.00 ±0.03 
0.40 ±0.04 

1430 ± 30 



d 
d 
d 



ms 
ms 
ms 
ms 
AU 



-1 
-1 
1 



Mo 
R0 
(cgs) 

Pq 
Mj 

Rj 
(cgs) 

Pi 

K 
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Table 2: Stellar parameters of WASP-22 from our 
spectroscopic analysis. 

Parameter Value 



^eff (K) 


6000 ± 100 


log 5 


4.5 ± 0.2 


6 (kms-1) 


1.2 ± 0.1 


vsini (kms~"^) 


3.5 ± 0.6 


[Fe/H] 


-0.05 ± 0.08 


[Mg/H] 


+0.06 ± 0.04 


[Si/H] 


+0.17 ± 0.11 


[Ca/H] 


+0.10 ± 0.11 


[Sc/H] 


+0.12 ± 0.12 


[Ti/H] 


+0.10 ± 0.07 


[V/H] 


+0.07 ± 0.09 


[Cr/H] 


0.00 ± 0.08 


[Mn/H] 


+0.02 ± 0.06 


[Co/H] 


+0.09 ± 0.13 


[Ni/H] 


+0.09 ± 0.06 


log N(Li) 


2.23 ± 0.08 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.2. 
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